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A  METHOD  OF  CALCULATING  THE  BIOLOGICAL  SHIELD  OF  LINEAR  PROTON 
ACCELERATORS  WITH  STRONG  FOCUSING 


V.  N.  Lebedev,  V.  V.  Mal'kov,  B.  S.  Sychev 

Formulas  for  calculating  the  biological  shield  of  linear  proton 
accelerators  with  an  energy  of  up  to  1  heV  [sic]  are  presented  in 
this  work.  A  plan  for  calculating  a  shield  is  developed.  A  graphical 
explanation  which  significantly  simplifies  the  estimate  calculations 
of  the  shield  is  given  for  the  basic  mathematical  formulas.  An 
example  of  calculating  the  shield  of  a  high-current  linear  proton 
accelerator  with  an  energy  of  800  meV  [sic]  is  shown.  The  presented 
method  allows  estimate  calculations  of  a  shield  to  be  conducted 
effectively  with  alternative  designing  of  the  accelerator's  structure. 

The  formalism  of  calculating  a  shield  differentiates  the  uni¬ 
formly  distributed  losses  of  a  beam  and  the  localized  impact  points. 

The  actual  compound  distribution  of  the  protons  which,  for  one  reason 
or  another,  emerged  from  acceleration  along  the  length  of  the 
accelerator  may  be  presented  in  the  form  of  a  specific  combination 
of  evenly  distributed  and  local  losses. 

Calculating  a  shield  in  the  case  of 
evenly  distributed  losses  of  a  beam. 

Taking  only  the  uniform  losses  of  a  beam  of  accelerated  particles 
into  consideration,  the  linear  accelerator  with  strong  focusing  may 

be  represented  by  an  infinite  anistropic  line  source.  The  line 

1 


density  of  secondary  penetrating  rad iat i or.-neutrons-of  that  source  is 
determined  by  the  expression : 


F(E) 


r^B  ev¬ 


en 


Here  Ip  is  the  proton  flux,  proton/s;  L  is  the  accelerator's  length; 

is  the  coefficient  which  character izes  the  value  of  uniformly 
distributed  losses;  B(Ep)  is  the  output  of  neutrons  per  proton  which 
emerged  from  the  acceleration  as  a  function  of  the  protons '  energy. 


The  output  of  neutrons  B(EP)  is  the  sum  of  evaporative  B^(E)  and 
cascading  Bj<(E)  neutrons  per  proton.  The  functions  B^(E)  and  B^(E) 
(Fig.  1)  are  plotted  on  the  basis  of  [1-3]. 


Fig.  1 .  The  dependence  of  the 
evaporation  and  cascading  neutrons' 
output  on  the  protons'  energy. 

KEY:  (1)  B(E)  neutron/proton;  (2) 

Proton  energy,  in  GeV. 


The  angular  distribution  of  evaporative  neutrons  (0.025  cV  <  E< 
<20  meV)  is  customarily  considered  isotropic.  In  accordance  with  thi 
we  write  the  evaporative  neutrons'  flux  depression  in  the  form  [4] 


fc\  Ji  ¬ 


lt  I 


(2) 


A  function  of  the  integral  secant  of  scci  (d/;\)  is  tabulated 
in  [4-‘3]  and  presented  in  Fig.  2. 


Fig.  2.  A  function  of  the  integral 
secant  of  P  (d , X  )  - 


In  the  formula  (2)  r  ii  the  distance  from  the  axis  of  a  proton 
beam  to  a  point  under  consideration  beyond  the  shield;  d  is  the 
shield's  thickness;  A  is  the  evaporative  neutrons'  relaxation  length 


The  evaporative  neutrons'  relaxation  lengths  are  not  determined 
experimentally.  Based  on  the  fact  that  the  evaporative  neutrons' 
spectrum  is  close  to  the  fission  spectrum,  we  used  the  relaxation 
length  for  neutrons  of  the  fission  spectrum ^  which  is  known  for 


The  fission  spectrum’s  A  ^  values  are 


many  materials  [6,  7]  as  A 
presented  in  Table  1  for  materials  which  found  application  in  the 
shield  of  proton  accelerators . 


We  write  the  law  of  flux  density  depress  i  on  of  cascading  neutron 
which  are  generated  by  an  anisotropic  liner  source  in  the  following 
manner 

-  jo.  (3) 

r 


Table  1.  Relaxation  lengths  for  neutrons  of  a  fission  spectrum  A  „ 
and  the  build-up  factors  of  B  and  B3.  & 
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KEY:  (1)  Material  in  the  shield;  (2)  Weight  by  volume,  kg/m3;  (3) 
Hydrogen  content,  kg/m3;  (4)  cm;  (5)  Build-up  factors;  (6)  Earth; 
(7)  Concrete;  (8)  Common;  (9)  Hematite;  (10)  In  steel  scrap. 


Table  2.  The  relaxation  lengths  of  cascading  neutrons  A  ,  in  cm. 
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(6)  Concrete;  (7)  Common;  (8)  Hematite;  (9)  In  steel  scrap. 


Formula  (3)  considers  the  angular  distribution  of  cascading 
neutrons  ((O)  [8],  the  analytical  expression  of  which  is  presented 
i n  the  form  [11]. 


2«f  (6)  =  0,2  f-  0.^3 cos ‘i  I, SI  cos* 9;  0<9  <— ; 

2 

— i-  9  it. 


('D 


2rf(0)~0.>  f  0,15  cos  9; 


The  expressi  on  which  stands  below  the  integral  sign  in  formula 
(3),  is  inter, rated  numerically.  The  results  of  the  computations  are 
presented  in  F i g .  3-  We  used  the  cascading  neutrons'  relaxation 
lengths  (Table  2)  from  [1]  for  the  calculations. 

Fig.  3-  Integral  function  of 

-d.  l-.ee  O/.T 

(ill  «  1  H  :i1.  d.  M  J  ' 


Calculating  a  shield  at  localized 
impact  points  of  a  beam 

The  localized  impact  points  of  a  beam  are  formed  in  areas  of  the 
placement  of  beam  gates,  Faraday  cylinders,  and  other  devices  which 
completely  or  partially  absorb  beams.  These  devices,  by  virtue  of 
the  small  dimensions  in  comparison  with  the  distance  to  the  point 
under  consideration  beyond  the  shield,  may  be  considered  point 
sources  of  neutrons  which  are  anistropic  for  cascading,  and  isotropic 
for  evaporative  neutrons. 
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The  flux  depression  of  evaporative  neutrons  from  the  source  is 
described  by  the  expression: 


(Ep,  r)  — 


(£')rti  -J>>( 

4rr-  C 


(5) 


The  flux  depression  of  high-energy  neutrons  (K  >  20  moV)  for 
these  same  conditions  may  be  presented  in  the  form 


<\(EP.  *. 


!p-BK(E-)r(,  -d> 

r)  =  ±—. - f{V  e 


(6) 


In  formulas  (5)  and  (6)  r  is  the  distance  from  the  source  to  the 
point  under  consideration  beyond  the  shield;  d  is  the  shield's  thick¬ 
ness;  f ( 0 )  is  the  angular  distribution  function  of  cascading  neutrons 
determined  by  formula  ( H ) ;  n2  is  the  coefficient  which  characterizes 
the  value  of  a  beam's  localized  impact  points.  The  remaining 
designations  are  the  same. 

A  scheme  for  calculating  a  shield 

s 

1.  Values  of  evenly  distributed  and  localized  losses  of  an 
accelerator's  beam  as  well  as  the  arrangement  of  localized  losses 
along  the  accelerator's  length  arc  determined. 

2.  The  entire  length  of  the  linear  accelerator  is  divided  into 
sections.  The  larger  these  sections,  the  more  exact  the  calcula¬ 
tion  of  a  shield.  It  is  quite  sufficient  to  select  5-6  sections  for 
e  s  t i ma  t  e  d  calculations. 

3.  The  densities  of  the  line  sources  (1)  and  (1)  are 
calculated  for  each  section  according  to  formula  (1)  and  with  the 
help  of  the  graphs,  (roe  Fig.  1).  It.  is  assumed  that  the  density  of 
the  source  is  constant,  along  each  sect. ion's  length. 
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I).  For  each  value  of  the  line  density  Fj  (1),  the  t  h  i  eknerr. 
the  shield  is  cal  eu  la  ted  accord  in;',  to  formula:;  (?)  and  (3)  b  mod  on 
the  ration: 

<t\i(r •  1',  dl  (l  •  If ,p)  --  (I>,  i?:  (7  ) 

<Mh  E.  <!)•(!  r  Bj)  —  ^n-j.y  (8) 

Here  Bp  in.  the  build-up  factor  of  intermediate  newt  re. nn. ;  Hv  is  the 
build-up  factor  of  decelerated  groups  of  radiations  (E  <  ?C  meV); 

<j>_  in  the  maximum  permissible  flux  densities  of  neutrons  beyond 

*  ■*  I 

the  s h i e  1  d  . 


The  values  of  B^p  and  B3  obtained  on  the  basis  of  [1,  9,  and  10] 
arc  presented  in  Table  1  and  the  recommended  values  of  On-y  -  in 
Tabic  3- 


Tabic  3-  The  recommended  values  of  the  max imum-permi 
fluxes. 
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KEY:  (1)  Group  of  neutrons;  (?)  Fast;  (3)  Very  fast  ;  ( A  1  furor  :'ar?  ; 
(!>)  Energy  of  the  rad  i  at  ion ;  (6 )  The  rated  flux  density,  n-*n:  ronr  /  cm 
for  places  .  .  .  ;  (6a)  which  are  basic  for  constant,  opera;;  us  with  ’la 
source  of  radiation  (catalyst.  A  —  1  )  ;  (bh)  which  are  rem  1  -re  *•  v  1  ■  b  . 
(catalyst  A-?);  (6c)  which  are  unserviceable  (catalyst.  A-,"’);  (od) 

which  are  adjacent  with  a  constant,  stay  of  pursonnol  (catalyst  B). 
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The  shield':-,  thickness  d  is  found  with  the  assigned  values  of  r 

from  relationships  (?.)  and  (  1 )  and  the  larger  of  d ,,  and  d  is 

*1  K 

taken  as  the  actual  value  of  d  in  each  section. 

5.  An  envelope  (curve)  which  is  also  an  estimated  outline  of  th 
shield  from  n  beam's  evenly  distributed  losses  is  plotted  according 
to  values  of  the  shield's  thickness  d  calculated  for  each  section  of 
the  accelerator. 

6.  The  shield  is  calculated  by  formulas  ( 'j )  and  (6)  in  areas,  of 
the  beam's  localized  impact  points;  a  final  outline  of  the  shield 
from  the  evely  distributed  and  localized  sources  is  plotted. 

An  example  of  calculating  an  800-meV  accelerator's  radiation  shield 

Initial  data.  Intensity  Ip  =  10iA|  proton/s;  energy  of  injection 
Ei  =  0.7  meV;  maximum  energy  of  the  protons  Ep  =  800  meV. 

The  acceleration  system  has  two  accelerating  structures.  The 
first  accelerates  the  protons  from  an  energy  of  0.7  meV  to  POO  meV, 
the  second  from  200  meV  to  800  meV.  The  lengths  of  the  accelerating 
structures  are  120  and  180  m  respectively  (see  Fig.  il ) . 


Fig.  H .  To  calculate  a  shield  with  evenly  dis¬ 
tributed  losses:  1  -  The  dependence  of  the 
protons.'  energy  on  the  length  of  the  aeeelerator; 
2  -  Calculated  outline  of  the  shield. 

KEY:  (1)  Sections;  (?)  Length  of  the  sections. 
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losses  n 
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Calculation.  1.  Wo  Laki-  the  value  of  the  evcly  distributed 
a  b  '  =  1  /•  ■ 


2.  We  divide  the  entire  length  of  the  accelerator  into  9 
sections.  The  length  of  each  section  is  60  n.  The  boundary  energy 
values  equal  0.7,  100,  200,  400,  600,  800  meV  (see  Fig.  4). 

3.  We  determine  the  output  of  evaporative  and  cascading  ncutrc-i 
for  each  of  the  accelerator's  section  by  the  graphs  (sec  Fig.  1): 


E  =r 

20  Mjti 

n„- 

io--. 

Bk  --<  I  a  -  ‘ . 

E  - 

100  M3U 

B,r  ' 

1.4  •  10-’. 

14  „  —  2,5  •  1 0  ’ 

E  =- 

200  .i i.Ki 

6.3  •  10 -l. 

Bk  --  2  •  1 0  - 1 ; 

E  — 

400  MJS 

1.7, 

Bk~  8,5-  10-’; 

E  = 

600  .1 no 

B„  = 

3.2. 

BK=r  1.6; 

E  — 
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5, 

2.5. 

4.  We  calculate  the  linar  source  densities  for  each  section  from 
the  formula  ( 1 ) : 


E  =  20  meV 


Fh  ~  3,3  -  !0‘  proton /m  •  s  . 

300 


F  w  'S  0; 


E  =  100  M3B 


P  1014 -0.01  •  1 ,4-  !0~l  f  . 

Fh  = - -  -4,7.10s  pro  ton /nr 

300 
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5  4 ■  1 Qr>  proton/m-s 
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Fg 
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5.  We  select  the  <t>.. ,  value  on  the  surface  of  Earth  which 
equal;'.  20  neulrons/cm'-  •  s .  We  find  the  build-up  factors 4  value  for 
Earth  from  Table  1  :  14  =  1  .2;  IC  =  3-9-  We  determine  the  permi  ss  i  M 

fluxes  for  evaporative  and  case ad i ng  neutrons  on  the  surface  of 
Earth  by  formulas  (7)  and  (8): 
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1 


4*., 


*J-  -  17: 
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<f», 

.vj 


6.  We  ca  1  cu  late  the  shield's  thickness  d(.  ,md  d/.  for  each 
section  o  i  the  a  c  o  e  1  e  r  n  t  oi1  L  y  formulas  (  2  )  a  n  d  (  S )  w  i  1.  h  the  help  of 
F  i  g  •’>  •  ?  and  3.  Upon  calculating,  we  assign  a  value  for  r,  then  we 
conduct  repetitive  calculations  for  the  purpose-  of  obtaining  the 
results  r  =  a  +  d,  where  a  is  the  distance  from  the  source  to  the 
shield. 

E  =  20  ineV 


We  assign  the  value  r  =  2.5  +  2  =  4.5  m.  Wc  determine  d 1 . 
From  formula  (2)  we  obtain  the  relationship: 

sed(a)'^-310-' 


From  the  graph  (see  Fig.  2)  we  determine  d„/A  =  1.2. 

h 

Table  1  =  18  cm;  consequently,  d  =  18*1.2  =  22  cm. 

A  Kt 

We  correct  the  value  of  r: 

r  =  2,5  -f-  0,5  —  3,0  .u; 


From 


17-4tt-300 


seci 


uj  3,3-10 


-  -  2*10-‘. 


From  Fig.  2  d ,  /A  =  1.4;  dM  =  18*1.4  =  25  cm.  The  value  d  -  0, 
r‘  jx  *'  la 

because  Fp(  =  0  when  the  pr'oton  energy  is  E  =  20  mcV. 

E  =  100  meV 


We  take  r  =  2.5  +  1.0 


5  m . 


We  determine  d 


11 


seci 


17  i-  :r.o 

■1,7 -10;  ip  • 


l,f.  1')--. 


10 


From  !•'  i  g .  2  d  =  3  ■  6 ;  d  =  18-  3.  6  =  86  cm. 
"  /L  17 

Wc  take  d  =  60  cm.  Wo  determine  d  . 

H  K, 

From  formula  (3)  we  obtain  the  relationship: 


\  U«)e  K  K 
0 

_ 

F|t(H)  8, 


■■  o 

<10  .  |  (<J,  < .  0;  <]'| 

6 


V-i.o 
■  H‘.  -  lC-  * 


2.1  10  6 


We  find  d  /A  =1.2  from  the  graph  (Fig.  3)-  When  the  energy 
is  E  =  100  meV  for  soil  \  =  69.1  cm.  Then  d..  =  69- 1- 1.2  =  70  cm. 

N  /N 

In  a  similar  manner,  we  determine  the  shield's  thicknesses  d v< 
and  d^  for  proton  energy  of  200,  J!00,  600,  and  800  meV.  The  values 
of  these  thicknesses  are  the  I'ol  lowing: 

E  =  200  mjj  d,,  =  95  cm,  dK  =  290  cm-, 

E  s=  400  M30  dM  --105  CM,  c)k“-  400  CM; 

E  —  600  M33  d.,  —  115  c.u,  dK— 450  c.u; 

E  — 800  M30  d„=120  cm,  dK— 460  cm. 


We  take  the  larger  values  of  d^  and  d^  as  the  actual  value  of 
the  shield's  thickness.  Those  are: 

E—  20  m:*3  d  =  25  c.u;  F  —  400  mjo  d  —  400  cm; 

E  =  100  mvs  d  —  70  c.u;  F.  —  600  M39  d  =  450  cm. 

E  —  200  .uaa  d~  290  c.u;  E  —  800  .ua<?  d  =  460  cm. 

We  plot  these  values  of  d  on  Fig.  >1  and  construct  the  estimated 

dependence  (2)  of  the  shield’s  thickness-  on  the  length  of  the 
accelerator  with  evenly  distributed  losses,  y  }  =  1%  of  the-  full 
intensity  of  Ip  ::  lO^1’  protons /s . 

1  3 
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